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ABSTRACT 

We report a detection of water in emission in the spectrum of the M2 
supergiant star /x Cep (M2Ia) observed by the Short Wavelength Spectrometer 
(SWS) aboard Infrared Space Observatory (ISO) and now released as the ISO 
Archives. The emission first appears in the 6 fim region (u 2 fundamental bands), 
and then in the 40 /iin region (pure rotation lines) despite the rather strong dust 
emission. The intensity ratios of the emission features are far from those of the 
optically thin gaseous emission. Instead, we could reproduce the major observed 
emission features by an optically thick water sphere of the inner radius about 
two stellar radii (2/2* « 13OOi? ), T cx = 1500K, and iV co i(H 2 0)= 3 x 10 20 cm" 2 . 
This model also accounts for the H 2 absorption bands in the near infrared 
(1.4, 1.9, and 2.7yum) as well. The detection of water in emission provides strong 
constraints on the nature of water in the early M supergiant star, and especially 
its origin in the outer atmosphere is confirmed against other models such as the 
large convective cell model. We finally confirm that the early M supergiant star 
is surrounded by a huge optically thick sphere of the warm water vapor, which 
may be referred to as MOLsphere for simplicity. Thus, the outer atmosphere 
of M supergiant stars should have a complicated hierarchical and/or hybrid 
structure with at least three major constituents including the warm MOLsphere 
(T ~ 10 3 K) together with the previously known hot chromosphere (T R3 10 4 K) 
and cool expanding gas-dust envelope (T R3 10 2 K). 

Subject headings: infrared: stars - molecular processes - stars: individual (/i 
Cep, a Ori) - stars: late-type - stars: supergiants 



- 3- 

1. INTRODUCTION 

Presence of water in the stellar environment has been known since 1960's, especially by 
the pioneering space observations with the balloon-borne telescope Stratoscope II (Woolf, 
Schwarzschild, & Rose 1964). However, it is the recent ISO mission (Kessler et al. 1996) 
that the infrared spectra unobscured by Earth's atmosphere have finally been explored in 
full details. One of the important results of the highly successful ISO mission is that water 
exists everywhere in the universe and that the star is not an exception. In fact, water has 
been detected not only by absorption but also by emission in cool luminous variables such 
as W Hya (Barlow et al. 1996; Neuffeld et al. 1996), NML Cyg (Justtanont et al. 1996), S 
Per (Tsuji et al. 1998), SW Vir (Tsuji, Aoki, & Ohnaka 1999), R Cas (Truong-Bach et al. 
1999), o Cet (Yamamura, de Jong, & Cami 1999), VY CMa (Neuffeld et al. 1999) etc. The 
exact nature of these emission features, however, is by no means clear yet. 

On the other hand, water in the early M (super)giant stars has been paid little 
attention until recently, since water has not been detected in these objects nor expected 
to be abundant in their atmospheres. Nevertheless, water in absorption has finally been 
detected by ISO in the early M (super)giants (Tsuji et al. 1997, 1998) as well as by 
the ground based mid-infrared observation in the early M supergiant stars (Jennings & 
Sada 1998). More correctly, however, water in the early M supergiants had in fact been 
discovered in the spectra taken by Stratoscope II (Woolf et al. 1964; Danielson, Woolf, 
& Gaustad 1965) more than 30 years before the ISO mission, although this discovery has 
been overlooked until recently in favor of the opposing interpretation. We have shown that 
the Stratoscope observers had correctly identified water in the early M supergiant stars but 
that the water detected by them is not originating in the photosphere (Tsuji 2000). 

To clarify the nature of water found in the early M supergiant stars, ISO should provide 
the best possibility. In fact, thanks to the recently released ISO Archives, we found water 
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in emission in \i Cep for the first time in such an early M supergiant star. This detection 
is somewhat unexpected, but this provides not only the definitive evidence for the presence 
of water in the outer atmosphere of the early M supergiant stars but also a new clue to 
understanding the nature of water in the stellar environment. 

2. OBSERVED EMISSION LINES OF WATER 

We retrieved the spectra of // Cep (M2Ia: observed on 1996 Dec. 18 UT) and a 
Ori (M2Iab: observed on 1997 Oct. 08 UT) from the ISO Archives. These spectra were 
observed with the SWS (de Graauw et al. 1996) by its highest resolution grating mode 
(resolution R fh 1600 or about 188 km s -1 ) in the region between 2.4 and 45 /mi. The 
region of the H 2 V\ and v% fundamentals around 2.7 /im is somewhat complicated due 
to the photospheric CO and OH bands, but the ISO data of both a Ori and \i Cep show 
additional absorption features not predicted by the photospheric models. At least, part of 
these features should be the same origin as the 1.4 and 1.9 \xm absorption bands, namely 
H2O of the non-photospheric origin (Tsuji 2000). 

EDITOR: PLACE FIGURE HERE. 

The region of the H 2 u 2 fundamentals (the band origin at 6.3 /zm) is shown in Fig. 1. 
First, the spectrum of a Ori shows only weak absorption features due to H 2 (indicated by 
the dotted lines in Fig. la), as can be confirmed by the comparison with the H 2 data based 
on HITEMP (Rothman 1997) shown in Fig.lc, but this detection is the best evidence for 
water in a Ori at present. The spectrum of \i Cep shown in Fig. lb is quite different with 
complicated spectral features which show upturn at about 6.4 /im in contrast to a smooth 
decline typical of the Rayleigh- Jeans spectrum as in a Ori. This upturn is explained by the 
onset of emission due to the P branch of the H 2 z/ 2 fundamentals. Thus, the H 2 bands 
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in absorption at 1.4, 1.9 and 2.7 /iin turn to emission at 6 /tin, and this can be regarded as 
definite evidence for the presence of water in the extended outer atmosphere of /i Cep. Also, 
some emission features of the observed spectrum are well identified with the fine structures 
of H2O bands by the comparison with the H2O data (Fig.lc). The intensity ratios of the 
emission features do not agree with those of the optically thin emission represented by 
the absorption cross-section of H 2 at T — 1500K in Fig.lc, and the situation cannot be 
improved by changing the temperature. 

The spectrum of \i Cep has once been observed by the NASA Airborne Infrared 
Observatory (Russell, Soifer, & Forrest 1975), and non-grey behavior from 5 to 8 /im has 
been noted. The observers preferred to assume excess absorption in the 4-5 /im rather than 
excess flux in the 5-8 /im. Another interpretation was to attribute the observed 5-8 /im 
feature as due to emission of H 2 z/ 2 fundamentals, and this was also proposing that water, 
besides dust, should be an important constituent of the outer atmosphere of /i Cep (Tsuji 
1978). There has been no means by which to test such a proposition at that time, but this 
is now fully confirmed by the ISO SWS which clearly resolved the fine structure of the H2O 
emission bands. 

We already know that H2O pure rotation lines were observed as absorption in a Ori by 
the ground-based high resolution spectroscopy in the 12 /mi region (Jennings & Sada 1998). 
However, the spectral resolution of the ISO SWS is not high enough to detect such faint 
lines. In the 40 /im region of a Ori, some emission-like features can be seen (Fig. 2a), but 
they cannot be identified with the known sources, except for an emission feature at 40.5 /im 
which roughly corresponds to H2O lines as can be known by the comparison with the H2O 
data (Fig.2c). However, it is still possible that the observed features are due to the mixture 
of absorption (including photospheric lines) and emission unresolved by the low resolution. 
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More distinct emission features appear in \i Cep (Fig. 2b) and most of them can be 
identified with the H 2 pure rotation lines by the comparison with the H 2 data (Fig. 2c). 
The intensity ratios of the emission features again suggest that the emission should be 
originating from an optically thick source. Such emission features are relatively well seen 
in the 40 region, but not so clear in the other mid-infrared region. One reason for this 
may be the presence of dust emission. To illustrate this, the entire spectrum of \x Cep 
observed by the ISO SWS is shown in Fig.3a by the red dots. It is clear that dust emission 
dominates the region longward of 8 //m. 

EDITOR: PLACE FIGURE | HERE. 

3. ORIGIN OF THE WATER EMISSION 

We showed that a hypothetical absorption slab of H2O with T ex = 1500K and 
Aral = 3.0 x 10 20 cm -2 above the photosphere explains the 1.4 and 1.9 fim bands observed 
in the Stratoscope spectrum of fi Cep (Tsuji 2000). The presence of emission at the longer 
wavelength in the spectrum of /1 Cep confirms that such a model of the non-photospheric 
origin of the water spectrum is basically correct, but we must now introduce geometrical 
extension of the water gas to explain the formation of emission features. For this purpose, 
we replace the simple absorption slab by a spherically symmetric isothermal gaseous 
envelope of water, whose inner and outer radii are Ri n and R ou t, respectively. The density 
distribution in this envelope is assumed to be in hydrostatic equilibrium under the gravity 
field of the central star, and hence major contribution to the spectrum comes from the 
relatively dense region near Ri n . For this reason, Ri n is more important than R ou t- 

We already suggested that the emission features may be originating from an optically 
thick source (Sect. 2). In fact, the column density (A co i = 3.0 x 10 20 cm -2 ) and the 
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absorption cross-section of H 2 typically 1CT 17 cm 2 mol -1 at T — 1500K (Figs, lc and 
2c) suggest that the optical depth of the strong lines should be as large as 10 3 . Thus, 
the water envelope is optically quite thick in the lines, and the warm molecular envelope 
can be referred to as molecular sphere, or MOLsphere for simplicity. Then, we must 
solve transfer equation in this MOLsphere and, since we are to interpret line spectra, the 
photospheric spectrum including the full details of the absorption lines should exactly 
be considered as the boundary condition. For this purpose, we apply the photospheric 
spectrum characterized by M* = 15M Q ,R* = 650 R®, L* « 6 x 1O 4 L and T cS « 3600 K 
(Fig.l of Tsuji 2000). The computation has been done with the spectral resolution of about 
50000 by the use of the H 2 line database HITEMP (Rothman 1997). After a few trials, we 
found that R in 2R* 13OOi? results in a spectrum that agrees with the observed one 
reasonably well. We also assumed R out ~ 4/?*, but this is of minor importance by the reason 
outlined above. The resulting emergent spectrum is shown by the green line in Fig. 3a. 

In the spectral region shortward of 5 fim, the MOLsphere contributes only to 
absorption, and the upper boundary of the spectrum shown by the green line in Fig. 3a 
represents the photospheric continuum, below which CO, OH, and SiO bands of the 
photospheric origin can be seen together with the H 2 absorption produced anew in the 
MOLsphere (molecules other than H 2 are not considered in the MOLsphere). At the 
shorter and longer wavelength sides of 6.3 fim, the R and P branches, respectively, of the 
H 2 z/ 2 fundamentals are seen in emission above the photospheric continuum (which is 
not explicitly shown but can be extrapolated from the near infrared continuum). This 
emission of water due to the MOLsphere remains prominent throughout the mid infrared 
region because of the stronger H 2 pure rotation lines. The MOLsphere contributes only 
to emission in the mid-infrared region and the lower boundary of the spectrum shown by 
the green line represents the photospheric continuum, below which OH pure rotation lines 
of the photospheric origin can be seen. But most of these OH lines are weaker or some OH 
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lines are missing as compared with those in the photospheic spectrum because of the filling 
in by the overlapping H 2 emission lines from the MOLsphere. 

For comparing the predicted spectrum from the MOLsphere discussed above with 
observations, we convolve it with the slit function of SWS which is assumed to be the 
Gaussian of FWHM = 187.5 km s _1 and the resulting low resolution spectrum is shown by 
the blue line in Fig. 3a. The result can directly be compared with the observed spectrum in 
the region shortward of 8 fim where dust emission is negligible, and is shifted vertically to 
fit the observed spectrum of \x Cep shown by the red dots. A close-up of the region around 
6.3 ^m is shown in Fig. 3b, where the red dots represent the observed spectrum and the blue 
line the predicted spectrum from the MOLsphere, as in Fig. 3a. The agreements between 
the observed and predicted spectra appear to be fine if not perfect. 

For interpreting the spectrum lonward of 8 /im, effect of dust emission on our predicted 
spectrum should be taken into account. Since dust envelope of /j, Cep is optically thin (e.g. 
Tsuji 1978), we simply add the dust emission to the spectrum by the MOLsphere. The 
region near 40^m is shown by an expanded scale in Fig.3c, where the red dots represent 
the observed spectrum and the black line the predicted spectrum which is the sum of the 
spectrum from the MOLsphere (the blue line in Fig. 3a) and dust emission represented by 
the Rayleigh- Jeans spectrum normalized to 162 Jy at 40 yum. The agreements between the 
observed and predicted spectra appear to be reasonable. This result also shows that the 
mid-infrared flux of apparently dusty stars such as fi Cep consists not only of dust emission 
but also of gaseous emission of water. The contribution of water to the thermal balance in 
the outer atmosphere as a whole is not very large (only 10 % or so in the case of [i Cep), 
yet water may be the dominant coolant in the inner part of the outer atmosphere. 

In the above analysis, we apply T ex and iV co i for H 2 in the MOLsphere determined 
by the previous analysis of the Stratoscope spectra (Tsuji 2000). We confirmed that 
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these parameters explain the ISO spectra as well and no improvement is gained if these 
parameters are changed. Certainly, these parameters are better determined from absorption 
rather than from emission spectra. Also, the near infrared spectrum predicted from the 
MOLsphere agrees well with the Stratoscope spectrum, which was interpreted by the simpler 
slab model in our previous analysis. The new constraint imposed by the ISO spectra is 
essentially the inner radius of the MOLsphere, Ri n . In fact, the emission features are quite 
sensitive to Ri n , and the emission tends to be much stronger for R in larger than R in ps 2R* 
we have adopted. Further, the 6 fim emission turns to absorption for the smaller values, e.g. 
R in ps lOOOi?© ps 1.5R*. The weak 6 fim absorption seen in a Ori (Fig.la) can be explained 
this way. Also, neither emission nor absorption will appear from the MOLsphere for some 
intermediate values of R in and the mid- infrared region of a Ori (Fig. 2a) may be such a 
case. Returning to fi Cep, the value of R in ps 2R* ps 1300/?© we have suggested should have 
important implications. Especially, this result suggests that water is absent in the region 
between r ps and 2i?*, and water of the MOLsphere may be formed anew somewhere at 
about R* above the photosphere rather than the region very close to the photosphere. This 
fact should be kept in mind in considering the origin of water in the MOLsphere. 

4. CONCLUDING REMARKS 

We found distinct emission of water at 6 fj,m and 40 jum regions in the ISO spectrum 
of n Cep in addition to the absorption bands at 1.4, 1.9 and 2.7 ^m in the near infrared. 
We showed that all these observed features can be explained consistently by a simple model 
of the optically thick molecular sphere with modest temperature (e.g. 1500K) at modest 
distance above the photosphere (e.g. about one stellar radius). What we have proposed, 
however, is not a model of the usual meaning but simply an intuitive interpretation on the 
observed data, since our "model" cannot yet be constructed based on the first principle. 
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In fact, the presence of the optically thick molecular envelope, which we referred to as the 
MOLsphere for simplicity, has never been predicted by any theory of stellar atmospheres 
so far as we are aware. Thus, despite the simplicity of the picture we propose, the physics 
behind the apparently simple phenomenon remains obscure and offers a challenging problem 
to the theory of stellar atmospheres. 

Recently, presence of water has been shown in several cool stars by ISO (Sect.l), 
but interpretation of the observed spectra is by no means clear. On the other hand, we 
could present a simple but consistent interpretation for the case of fi Cep and this may 
be due to fortunate situations as follow. First, n Cep is unique in that water could be 
observed through the near- to mid-infrared (including the Stratoscope data) and not only 
in absorption but also in emission. Second, \i Cep has the MOLsphere which has just 
developed to an adequate size. For comparison, a Ori may be the case in which the 
MOLsphere is not yet developed enough to give clear observed features, or emission and 
absorption (including the photospheric contributions) may be nearly canceled. On the 
contrary, the cooler supergiants may have well developed MOLsphere, the observed features 
from which may be heavily saturated and more difficult to interpret. Third, the spectrum of 
\i Cep shows almost no photospheric water band, unlike the case of the cooler objects which 
show strong photospheric water bands difficult to discriminate from those originating in the 
MOLsphere. Fourth, the MOLsphere of fj, Cep may be relatively free from complexities due 
to the dynamical effect such as by a large amplitude pulsation. For all these reasons, the 
spectra of \i Cep could be interpreted relatively well and this may provide the proto-type of 
the MOLsphere, which should most probably exist in the outer atmosphere of cool luminous 
stars in general. 

We finally confirmed that water exists not only in very cool luminous variables, but 
also in the early M supergiants and possibly in most M (super)giant stars. This fact implies 
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that the presence of water is not due to some secondary effect such as a large amplitude 
pulsation but may be due to a more intrinsic property of cool stars. For example, the 
MOLsphere may have some connection with the chromosphere which also exists in the outer 
atmosphere of the most cool stars, but its origin is by no means clear yet. By the opening 
of the infrared spectral region by ISO, we have now a means by which to explore the new 
component of the outer atmosphere which is rather warm (MOLsphere), in addition to the 
hot chromosphere and the cool expanding gas-dust envelope known so far. The problem 
is how to understand these multiple components of different thermal and kinematical 
properties in a unified model, and we now have a possibility to arrive at a more balanced 
view of the stellar outer atmosphere by considering all its major constituents. 

I thank W. Aoki, T. Tanabe, and I. Yamamura for their help and useful discussion on 
the ISO data. This research is supported by the Grant-in- Aids for Scientific Research of the 
Ministry of Education, Culture, and Science No. 11640227. 
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Fig. 1. — (a) Spectrum of a Ori by the ISO SWS (resolution R m 1600) and corrected for the 
interstellar reddening with A v = 0.5 mag (Lee 1970). (b) The same but for /x Cep corrected 
for the interstellar reddening with A v = 1.5 mag (Lee 1970). (c) Absorption cross-section 
(cm 2 molecule" 1 ) of H 2 at T = 1500K based on HITEMP (Rothman 1997). The black line 
is by high resolution (R 50000) and the white line by low resolution (R R3 1600). 

Fig. 2. — (a) Spectrum of a Ori. (b) Spectrum of \x Cep. (c) Absorption cross-section of 
water at T — 1500K. As for details, see the legend to Fig.l. 

Fig. 3. — (a) The full spectrum of /j Cep by the ISO SWS (corrected for the interstellar 
reddening) is shown by the red dots. The gaps around 4, 7, 12, 20, 29 /im are instrumental 
effect not yet corrected for. The predicted spectrum based on the MOLsphere (optically thick 
isothermal sphere of the water gas with T cx = 1500K, A" co i = 3.0 x 10 20 cm -2 , and located 
at .R* = 65O-R above the photosphere) computed with the photospheric spectrum as the 
boundary condition is shown by the green line (resolution R m 50000) and by the blue line 
(R « 1600), after being shifted vertically to fit the observed spectrum. For comparison with 
observation, the effect of dust emission should be added for A > 8/iin. The molecular bands 
indicated are originating in the photosphere except for H2O formed anew in the MOLsphere. 
(b) A close-up of Fig. 3a around H 2 v 2 fundamentals. The observed spectrum is shown 
by the red dots (identical with the ones in Fig. 3a) and the predicted one based on the 
MOLsphere by the blue line (identical with the one in Fig. 3a). (c) The region of the H 2 
pure rotation lines in an expanded scale. Now the dust emission is added to the predicted 
spectrum from the MOLsphere (the blue line in Fig. 3a), and the resulting spectrum is shown 
by the black line for comparison with the observed one (the red dots). 
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